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A polaron is an electron interacting with a polar crystal, which is able to form a bound state by 
using the distortions of the crystal induced by its own density of charge. In this paper we derive 
Pekar's famous continuous model for polarons (in which the crystal is replaced by a simple effective 
Coulomb self-attraction) by studying the macroscopic limit of the reduced Hartree-Fock theory of 
the crystal. The macroscopic density of the polaron converges to that of Pekar's nonlinear model, 
with a possibly anisotropic dielectric matrix. The polaron also exhibits fast microscopic oscillations 
which contribute to the energy at the same order, but whose characteristic length is small compared 
to the scale of the polaron. These oscillations are described by a simple periodic eigenvalue equation. 
Our approach also covers multi-polarons composed of several electrons, repelling each other by 
Coulomb forces. 

© 2011 by the authors. This paper may be reproduced, in its entirety, for non-commercial purposes. 
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Introduction 

In vacuum, N electrons cannot form a bound state, because of their Coulomb repulsion and 
the dispersive nature of their kinetic energy. Confining the electrons in a given region of 
space is only possible by applying an external field. In atoms and molecules, this field is 
the electrostatic attraction of the (positively charged) nuclei. The situation is completely 
different when the electrons are placed in a polarizable medium like a dielectric crystal. 
There, the N electrons induce a lattice distortion by repelling a bit the negative ions of the 
crystal and attracting the positive ones. When the resulting polarization is strong enough, 
the electrons are able to overcome their Coulomb repulsion and to form a bound state, called 
an N -polaron. 

Polarons have been widely studied in the physics literature. The main difficulty is to ad- 
equately describe the behavior of the polar crystal and of its distortions. The simplest model 
was obtained by Pekar [39l [40] who proposed to approximate the crystal by a continuous 
polarizable medium, described only by its static and high frequency dielectric constants. For 
the polaron, composed of only one electron, one gets the energy functional [H [5] 

£? M bP,D] = l [ m{x)fdx-f \i,(x)\ 2 v(x)dx+ 1 f \D(x)\ 2 dx. (0.1) 

Z JR3 Jm. 3 S7T^(£m] -tJJ»3 

Here ip is the wave function of the electron, D = —Vv is the displacement field of the 
continuous medium, £m is its static dielectric constant, in units such that the high frequency 
dielectric constant is Eq = 1 and such that the charge and the mass of the electron are 
normalized to e = 1 and m e - = h 2 . For simplicity we neglect the spin of the electron in the 
whole paper. Minimizing the above functional with respect to D at fixed ifi, we obtain the 
following energy functional for the polaron alone in the continuous polarizable medium 



CM - \ I |v*WI>* + (SMtini f [ (o.2) 

2J K 3 2 Jr3Jrz \x-y\ 



The last nonlinear term in (|0.2I) is an effective Coulomb self-interaction. When £m > 1, 
this nonlinear term is attractive and it has been shown by Lieb [27] that the Pekar energy 
functional (10.21) admits a unique minimizer up to translations (under the normalization 
constraint J R3 |-0| 2 = 1). This ground state is radial and solves the corresponding Eulcr- 
Lagrange equation 

(-y + ((^m)- 1 - * M" 1 ) $ = E1>. (0.3) 

This nonlinear equation is ubiquitous in Physics and it is sometimes also called the 
"Choquard" or "Schrodingcr-Newton" equation. 

Pekar's theory can be easily generalized to the case of N electrons, as was first suggested 
for N = 2 by Pekar and Tomasevich [41]. Taking into account the electrostatic repulsion 
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between the electrons, one arrives at the following many-body energy functional: 




The many-body wave function ^ must be antisymmetric with respect to exchanges of the 
variables x\, ...,xn, due to the fermionic nature of the electrons. Also, is the total density 
of the N electrons, defined by 



There is now a competition between the many-body electronic repulsion and the nonlinear 
attraction due to the polarizable medium. It has been shown recently by one of us [25] that 
when £m is sufficiently large (depending on TV), the many-body Pekar functional (|0.4p admits 
at least one minimizer, hence infinitely many by translation invariance. On the other hand, 
it can be deduced from the results of Frank, Lieb, Seiringer and Thomas jT5J Q3] that when 
£m < 1 + a (with a independent of N), £f M has no ground state for N > 2. 

Pekar's functional is not the only one used by physicists to describe (iV-)polarons in- 
teracting with a continuous medium. In [16j Q~Z] > H. Frohlich has proposed to replace the 
classical polarization field D of (jO.ip , by a quantized (phonon) field with which the electrons 
interact [2]. Frohlich's model has been mathematically studied in several works. In partic- 
ular, Donsker and Varadhan [12] and then, with a different approach, Lieb and Thomas 
[32] . have proved that Pekar's polaron can be recovered from the strong coupling limit of 
Frohlich's model. This was later extended to bi-polarons by Miyao and Spohn in [36] ■ For 
other recent works on Frohlich's and Pekar's theories, see for instance [37 ] [T9 ] rT5 ] IM ] IS]. 

Both models assume that the medium in which the particle evolve is continuous. In a 
crystal, this can only be valid when the size of the electronic system is much bigger than 
the typical lattice length, that is, the diameter of the unit cell. One then speaks of large 
polarons. For smaller polarons, this approximation is not good enough and the electrons 
start to see the detailed structure of the crystal. 

The purpose of this paper is twofold. First, based on previous works by the first author 
with Cances and Deleurence [8] [9] [10] , we write a simple mean- field model describing small 
multi-polarons in a quantum crystal. Second, we show that in a macroscopic limit where 
the system lives on a scale much larger than the size of the lattice cell, we again recover the 
Pekar(-Tomasevich) theory. We now explain the main lines of our approach, before turning 
to a more detailed presentation of our results in the next section. 

Our crystal is assumed to be extended over the whole space. At rest, it is composed 
of classical nuclei, described by an Jzf-periodic charge density /Zp er > 0, and of quantum 
electrons which are modelled by an Jzf -periodic density p pcr . The lattice Jzf is a discrete 
subgroup of M 3 whose fundamental domain T (the unit cell) is compact, for instance Jzf = Z 3 . 
The system is locally neutral in the sense that J r Pp Gr = J r Mp 0r - The unperturbed crystal 
induces an Jzf-periodic electrostatic potential Vp° or which solves Poisson's equation 



This electrostatic potential is felt by any other particle which is added to the system. 




(0.5) 
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When the additional particles are inserted, the nuclei and the electrons of the crystal 
can be displaced a little bit. This distortion is described by (local) perturbations Sp and 
Sp such that the nuclear and electronic densities become p = Pp el + Sp and p = /?p Cr + dp. 
The inserted particles then feel the electrostatic field (Sp — dp) * | | 1 induced by these 
displacements. For the single polaron, we are thus led to an energy functional of the form 

£[ip,5p,8n\ = \ [ \V^(x)\ 2 dx + I V£ ei {x)\^(x)\ 2 dx 

— ?^)——J^ j) 1^(^)1 dxdy + FciysiSpjSp] (0.6) 
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where J- C ry S [Sp, Sp] is the energy cost to perturb the crystal by moving the nuclei of Sp and 
the electrons of Sp. Note that the periodic density p® cr — p® CI is locally neutral, and that 
the displacement densities Sp and Sp should satisfy 

Sp(x) dx — I 8p(x) dx = 0, 

JR 3 

at least formally (see Remark 1 1 . 1 1 below) . As expected, the polaron effectively sees an elec- 
trostatic potential induced by a field of dipoles. 

Our two subsystems (the polaron and the crystal) are uncorrelated in this simplified 
theory, hence it is possible to completely eliminate the crystalline degrees of freedom, by 
minimizing over Sp and Sp for any fixed state tp of the polaron. This leads to an effective 
nonlinear functional for the electron alone, of the form 

£*sbP} = l [ \Vi>(x)\ 2 dx+ ( V° er (xM(x)\ 2 dx + F ciys [\ilj\ 2 ], (0.7) 

z JR 3 JR 3 

where the nonlinear effective energy F crys is defined by 

iWH 2 ] = . -f n [I I ^ ~^ A l C ^ )l ^ C2/)|2 • (0.8) 

The case of the iV-polaron is obviously similar, leading to the effective nonlinear many-body 
functional 




\Xk — Xf\ 
Kk<i<N ' K 



+ / V° ei {x)py(x)dx + F ays [py], (0.9) 

where -Fcrys is the same nonlinearity as in (|0.8[) . 

Of course, the main difficulty in this context is to find an energy functional J^cryslSp, Sp] 
quantifying the cost to move nuclei and electrons in the crystal, which is both physically 
relevant and mathematically amenable. It turns out that there is a well-defined mean-field 
theory [8] for the electronic perturbation Sp for every fixed value of Sp, but that allowing the 
nuclei to move freely is too involved for the present mathematical technology. If the crystal 
is not globally stable with respect to the positions of the nuclei, the latter will want to relax 
to better positions, changing thereby the total energy per unit volume and rendering our 
above effective energy -Fc-ys infinite. 
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To face this problem we could add some stability conditions on the nuclear structure 
of our crystal, but this would complicate our exposition dramatically. Since in this paper 
we are more interested in the derivation of Pekar's polaron from a microscopic model than 
in proposing a quantitative theory, we will make the (very strong) assumption that the 
nuclei cannot move, <5/i = 0. As we will see, the distortion of the electronic Fermi sea of the 
crystal is in principle enough to bind polarons, although its effect is weaker than when the 
nuclear displacements are taken into account. As far as the derivation of Pekar's polarons 
is concerned, this simply means that we will only obtain the electronic contribution to the 
dielectric constant em- 

In Section 11.11 below, following [8] , we define an appropriate functional -F C rys obtained in 
a Hartree-Fock-type approximation, when only the electronic Fermi sea is allowed to move. 
In this introduction we assume that F crys is given to us without giving its precise expression, 
and we now discuss the derivation of Pekar's polarons in a macroscopic limit. Of course the 
precise form of -F crys is very important, since it is at the origin of the dielectric constant £m 
seen in Pekar's theory. 

We now explain our derivation of Pekar's energy functional in the case of only one electron 
(polaron), the argument being similar for the ./V-polaron. We have to let the polaron live 
on a much larger scale than the typical size of the lattice cell or, equivalently, to make the 
crystal live on a much smaller scale than that of the polaron. To this end, we introduce a 
small parameter < m< 1 which is interpreted as the ratio between the microscopic and 
the macroscopic lengths. At the macroscopic scale, the lattice becomes mJ*? and the energy 
functional of the polaron is now 

£mbP] = l [ \S7^{x)\ 2 dx + m- 1 [ V^ix/m^ix^dx + m^F^m^im-)] 2 ]. (0.10) 

* JR 3 JR 3 

The scaling is chosen to make all the terms of the energy contribute the same at the macro- 
scopic level. As for the periodic potential of the crystal, we have that V m :— m^ 1 Vp Cr (x/m) 
is the unique solution of Poisson's equation 

-AV m = m~ 3 (p° pcl (x/m) - ju£ er (x/m)) . 

The scaling of the density is chosen such as to keep constant the number of electrons 
and nuclei per unit cello There are similar arguments in favour of the chosen scaling 
™ 1 F cryB (m 3 \il}(m-)\ 2 ) for the nonlinear term. 

Changing variables ip = m 3 / 2 ijj(m-) we can express the same functional at the microscopic 
scale 

£ m $\ =m- 1 ^^jV^(x)| 2 dx + £ 3 y p ° cr (x)|^(x)| 2 rfx + ^ crys [|V;| 2 ]^ . (0.11) 

From the perspective of the crystal, a large polaron can therefore be obtained by inserting 
a particle whose mass m is very small, and which thus tends to be very spread out in space. 

For an isotropic crystal described by the nonlinear energy F C rys defined later, we prove 
in Theorem 11.11 below that in the limit m — > 0, any ground state of (|0.10j) behaves as follows 



(0.12) 



lr That is, we have f mF m 3 p£ cr (7 m ) = / r p° eI and j mV m 3 /x° cr ( /m) = / r n° e 
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up to a well-chosen translation of the system in space. Here ip^ is the unique ground state 
of Pekar's functional (|0.2[) . £m > 1 being the macroscopic dielectric constant of the crystal, 
which will be defined in Section 11.11 On the other hand, is an Jzf-periodic function, 
which converges to 1 uniformly as m — ¥ 0. It is defined by minimizing the functional 




with periodic boundary conditions on dT and under the constraint that L |uP^ r | 2 = 
where we recall that T is the unit cell of the lattice. Extended by periodicity over the whole 
of R 3 it solves 

(~ + Cw) < cr ( a; ) = E ™ r (°- 13 ) 

The precise behavior of uf^* as m — > (as well as the value of Ef") can be determined by 
usual perturbation theory, as we will explain later in Section [5] The corresponding energy 
of ip m tends to the sum of the energies of the two functions u^ r and ip^ M : 

Urn E m [An] = + £f M [^J (0.14) 

where E pcl is the limit of ra^ 1 E^ T when m — > and £f is Pekar's energy (|0.2[) . 

Let us emphasize that the limit studied in this paper is completely different from existing 
results on Frohlich's model [T^l GUI [35] . In Frohlich's theory, £m is a parameter to be chosen. 
The goal is to show that in the strong coupling limit the polaron tends to decouple from 
the quantized field, leading to Pekar's ground state with the given £m- In this paper the 
correlations between the polaron and the crystal are already neglected and our purpose is to 
derive Pekar's model in a macroscopic limit. Our derivation provides a certain value for the 
dielectric constant £m, in terms of the structure of the chosen microscopic quantum crystal. 

In this introduction we have explained the simplest situation of an isotropic crystal whose 
dielectric tensor £m is a constant. Below we consider the general case and, in the anisotropic 
case, we obtain a generalized Pekar functional in which £m is a 3 x 3 real symmetric matrix. 
The corresponding expression for £f M will be given in Section 11.21 below. Also, the results 
discussed here for the single polaron hold similarly for iV-polarons. 

The paper is organized as follows. In Section [T] we describe our small polaron model in 
detail. We then recall some facts about Pekar's theory and state our main results relating 
the latter to the reduced Hartree-Fock theory of quantum crystals. Their proofs rely on 
two main ingredients. First, in Section[5]we separate out the contribution of the microscopic 
oscillations, by using a simple energy decoupling argument and the properties of the periodic 
eigenvalue problem. The second, more involved, step is the detailed analysis of the limit of 
the perturbed crystal model of jS]. Section [3] proceeds with improving some results of [10] 
and applying them to the context of the polarons. Finally, we complete the proofs of our 
main results in Section [4] 

Acknowledgment. The research leading to these results has received funding from the Eu- 
ropean Research Council under the European Community's Seventh Framework Programme 
(FP7/2007-2013 Grant Agreement MNIQS No. 258023). 
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1. Main results 

In this section we introduce some preliminary tools which are necessary to properly define 
our model, and we state our main theorems. Proofs will be given in the next sections. 

1.1. Energy to perturb the interacting periodic Fermi sea 

In this section we quickly present the model that Cances, Deleurence and Lewin have in- 
troduced in [8j [9] to describe the distortions of the electrons in a quantum crystal, using 
a Hartree-Fock-type theory. We will define the effective nonlinear energy F crys felt by the 
polaron, and also quickly recall the definition of the macroscopic dielectric tensor em which 
will later appear in Pekar's functional. We provide some additional technical details later in 
Section [3TT1 

We fix an Jzf-periodic density of charge /ip er for the classical nuclei of the crystal, with 
a discrete subgroup of M 3 . It is enough for our purpose to assume that fip CI is a locally- finite 
non-negative measure, such that J r ^tp C1 . = Z G N, where r = R 3 /Jzf is the unit cell. 

In reduced Hartree-Fock theory [45] , the state of the electrons in the crystal is described 
by a one-particle density matrix, which is a self-adjoint operator 7 : L 2 (R 3 ) — > L 2 (R 3 ) such 
that < 7 < 1 (in the sense of operators). When no external field is applied to the system, 
the electrons arrange in a periodic configuration 7 = 7p er , which is a solution of the reduced 
Hartree-Fock equations 

'7^r = l(-ac,e P )(-A/2 + y p ° er ), 
< - A ^p°cr = 47r (P7°„ r " Mp er )> (1.1) 

Here pa denotes the density of the operator A which is formally given by pa{x) = A(x,x) 
when A is locally trace-class. Also, l(_ 00)Ep )(ff) denotes the spectral projector of H onto 
the interval (—00, ep). The real number ep in Ijl.ljl is called the Fermi level. It is also the 
Lagrange multiplier used to impose the constraint that the system must be locally neutral 
(third equation in (|1.1|) ). It should not be confused with the dielectric constant £m which will 
be defined later. Existence and uniqueness of solutions to the self-consistent equation (|1.1| 
was proved in [H][5] (see also [35] at positive temperature). 

By the Bloch-Floquet theory [32] , the spectrum of the Jzf-periodic Schrodinger operator 

is composed of bands. When there is a gap between the Zth and the (Z + l)st bands, the 
crystal is an insulator and e-p can be any arbitrary number in the gap. Like in [5], in the 
whole paper we will assume that the host crystal is an insulator. 

Assumption 1.1 (The host crystal is an insulator). 

The periodic Schrodinger operator H® cr has a gap between its Zth and (Z + l)st bands, and 
we fix any chemical potential sp in the corresponding gap. 

When the quantum crystal is submitted to an external field, the Fermi sea polarizes. The 
new density matrix 7 of the system now solves the nonlinear equation 

7 = l(-oo, eF ) ( - A/2 + y p ° cr + (p T _ 7 o M - v) * III" 1 ) + S, (1.2) 
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where v denotes the external density used to perturb the Fermi sea, and where the operator 
S (satisfying < S < 1) lives only at the Fermi level!. Existence of solutions to this equation 
was shown in |S]. In general there is no uniqueness of 7, but the density p 7 is itself unique. 
The perturbed state 7 is such that Q :— 7 — 7 p er is Hilbert-Schmidt and locally trace-class, 
and pq G L 2 (M. 3 ). In general, Q is not even a trace-class operator |10) . This motivates the 
introduction of a particular functional setting, some of its elements being recalled below (see 
also Section RTTj) . 

The method used in [8] to construct solutions is variational and it relies on an energy 
functional which we want to use in our polaron model. The idea is to define the energy cost 
to move the electrons from 7 p er to 7 as the (formal) difference between the (infinite) reduced 
Hartree-Fock energies of 7 and of 7 p er . Denoting by 



DU , g) .. = ff mm^.^f /wwfi dk (1 .3, 

JJrsxw \x-y\ J R 3 \k\ z 

the Coulomb interaction, the formal calculation is 

£ rHP M - £ rHF [7p°cr] " = " Q Tr(-A - e F ) 7 + \d{ Pi - Mper , p- t - ^ c 

\ Tr(-A - £ F )7per + \ D (PlU - Mpcr, P-y°„ - Mper) 

" = " Tr(ff° er - e F )Q + \d{ Pq , Pq ) (1.4) 

with again Q = 7 — 7p 0r . Note that Q satisfies — 7° < Q < 1 — 7p Cr , which is equivalent 
to Q 2 < Q++ - Q— with Q++ := (7^)^(7^ > and Q~ := 7 ° cr Q 7 ^r < 0. This 
allows to properly define the kinetic energy in (jl.4p as follows |S]: 

Tr (ff p ° er - e F )Q = Tr \H* cl . £f | 1 / 2 (Q++ - Q~)|flJ„ - £ f| 1/2 . (1.5) 

More generally, one can define the generalized trace as 

Tr Q = Tr Q ++ + Tr Q (1.6) 

when Q ++ and Q are trace-class. 

The relative energy (11.51) is > for every e F in the band gap. The total energy to go 
from 7 p er to 7 is defined as 

•Fcry S [Q] := TroOffp° cr - e F )Q + ^D(p Q ,p Q ) (1.7) 
and it is also non-negative since we have 

2 



«r 



Taking into account the term involving the chemical potential e F in the total energy allows 
to have a non-negative energy functional (the reference periodic Fermi sea now has energy 
zero). The precise value of e F in the gap does not matter at this stage. 



2 This means Ran(<5) C Ker( — A/2 + Vp Cr + (p 7 _ 7 o — l/)*|x| 1 — ep ). The operator 5 can safely be ignored 
by the reader, as in the macroscopic limit that we will later consider, we will always have 5 = 0. 



Derivation of Pekar's Polarons 



9 



When we submit the crystal to an external density v, the state of the Fermi sea is 
obtained by solving the following minimization probler 



F CIys [i/} = inf [D{v, p Q ) + J" crys [Q] 



n| 

^ (1-8) 



As shown in [5], for any v 6 L 1 (K 3 ) n L 2 (K 3 ), this minimization problem has at least one 
solution Q = 7 — 7° (in an adequate function space that we recall in Section with 
7 solving the nonlinear equation (|1.2p . The corresponding density pq is in L 2 (R 3 ) but in 
general it has long range oscillations which are not integrable at infinity [TU]. The total 
energy i^crysM of the crystal submitted to an external density v is the one which we will 
use later, with v = \ip\ 2 (polaron) or v = p^ (7V-polaron). 

Remark 1.1. Instead of working in a grand canonical formalism where £f has a fixed value, 
we could impose that no charge can be added to the crystal which, in this context, means 
Tro(Q) = (recall that pq is not necessarily integrable at infinity). It has been shown in 
(Lemma 5) that when the external density v is small enough, e.g. D(v,v) < rj, then we have 
automatically Tro(Q) = for every ep which is at a distance > C \Jrj of the edges of the gap. 
The value of F crys [v) is independent of sp in this range. 

In the macroscopic limit that we consider later, the polaron induces a charge defect in the 
crystal which, at the microscopic scale, is very spread out in space and we will automatically 
have that the Fermi sea stays neutral, Tro(Q) = 0. As for the derivation of Pekar's polaron, 
the precise value of Ef in the gap therefore does not matter. 



Remark 1.2. Using well-known ideas from Density Functional Theory \28f . one can express 
the whole problem only in terms of the perturbed density 5p — pq, as we have done in the 
introduction. The energy cost to perturb the Fermi sea by a density bp is defined as 

•^c'rysIM = n inf n -^crysIQ], 

-7pc,-<Q<l-7p„ 

PQ=Sp 

such that -FcrysM can also be expressed as 

F CIys [v] = inf [D{v, p Q ) + T a s [Sp 

When the crystal is submitted to the external density v, the corresponding total energy 
-FcrysM satisfies the following simple estimate. 

Lemma 1.1 (A uniform estimate on F crys ). 

We have for all v G L X (R 3 ) n L 6 / 5 (R 3 ) 

-^D(v,u)<F ciys [u}<0. (1.9) 

Proof. The upper bound is obtained by taking Q — as test state in (|1.8[) . For the lower 
bound we neglect the positive kinetic energy of Q, complete the square and use that D{-, ■) 
defines a scalar product. □ 



3 In the whole paper we use a sign convention for the external density v which is opposite to that of |8lll0|. 
In these works v was interpreted as a nuclear defect density, whereas in our case it will be that induced by 
our polaron, which are negatively charged particles. 
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The estimate (|1.9j) already has interesting physical consequences for the polaron. The 
fact that -Fcrys is non-positive means that the interaction of the electron with the crystal will 
always be attractive. On the other hand, the lower bound means that Pekar's interaction 
(with em = oo) is always a lower estimate to the total interaction with the crystal. 

However, we will see below that (jl.9p is not optimal in the macroscopic limit, correspond- 
ing to density z/s which are very spread-out in space. The correct macroscopic behavior was 
studied by Cances and Lewin in |10) . In particular, they have derived the (electronic) di- 
electric matrix em of the crystal, which is a 3 x 3 symmetric real matrix such that £m > 1 
in the sense of matrices. 

The dielectric matrix em can be expressed in terms of the Bloch transform of i?p Cr , see 
Eq. (36) in [TO]. The formula of £m is well known [I] [48] but it is not really important to us. 
More important is Theorem 3 of [10j . which says that £m can be obtained by a macroscopic 
excitation of the Fermi sea: If we take an external density of the form 

v m(%) — m3 v(mx) 

and call Q m one corresponding solution of (jl.2p . then the rescaled self-consistent potential 

W m := m- x {y - p Qm ) * \ ■ ^(x/m) (1.10) 

converges weakly to the unique solution W v of Poisson's equation 

- div(e M VW^) = Airv. (1.11) 

The matrix em is the one which will appear in Pekar's theory below. The main reason for 
this is the fact, not derived in |10] but proved in this paper, that the energy behaves as 

= 27T L |9(fc)|2 G4* - 1^) *■ ^ 

for any fixed is, see Theorem II .41 below. 



1.2. Pekar's polarons in an anisotropic continuous medium 

In this section we introduce the Pekar functionals which will describe the macroscopic behav- 
ior of our (iV-)polaron. Since we want to consider anisotropic crystals, we have to generalize 
the formulas quoted in introduction (|0.2p and (|0.4I) to the case of Em being a 3 x 3 symmet- 
ric matrix. In the macroscopic limit considered later, the matrix em > 1 will be that of the 
reduced Hartree-Fock crystal derived in [JO] and recalled in the previous section. 
We start with Pekar's generalized functional which we define as 

CM : = J / m 2 dx + F^\ 2 ] (1.13) 

with 

Alternatively we can define £f M and F^ M by considering the solution W p of Poisson's equation 

- V (e M VWp) = 4np. (1.15) 

We then have 

£lM-=\l \vM 2 dx + \( M 2 - M 2 * I • I -1 ) • (i-ie) 
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We define the corresponding ground state energy as 

The following deals with the existence of minimizers for the variational problem (|1.17j) . 

Theorem 1.1 (Ground states of Pekar's functional, N = 1, [27], 133] ) ■ 

We assume that em > 1 in the sense of symmetric matrices. 

(1) (Existence). The minimization problem (| 1 . 1 T|) admits at least one minimizer ipf M - H 
satisfies the nonlinear equation 

l^L + (w^r - ICI 2 * i • r 1 ) C = a C- (i-is) 

(2) (Convergence of minimizing sequences). All the minimizing sequences (ipn) for 
Ef (1) are precompact in H 1 (M 3 ), up to a translation. That is, there exists a sequence 
(rfc) C R 3 and a minimizer ijjP of Ef (1) such that i/j nk (' ~ T k) ~ * ^ P M strongly in 

This result is a standard statement on which we will not elaborate. It can be proved by 
using well known techniques of nonlinear analysis, which are similar to what we do later 
in Section 13.41 When £m is proportional to the identity, Theorem 11.11 was proved first by 
Lieb in |27j . and in this special case the minimizer is unique (up to translations). Other 
functionals very similar to £f M have been considered by Lions in [33J for N = 1. To our 
knowledge, uniqueness for an anisotropic em is not known. 

Of course, the main reason why Ef (1) always has a minimizer is that the resulting 
potential is attractive at long distances. This is the case even when sm has 1 as an eigenvalue, 
as soon as £m 1 (we have k T eyik > \k\ 2 except on a set of measure zero). 

We now turn to Pekar's multi-polaron problem, whose energy functional in an anisotropic 
medium is 

J R 3« ^ j = 1 l<k<l<N \ Xk Xl \ 

+ F £ P >*] (1.19) 

where Ff is defined above in (|1.14[) and with p$ the density of the many-body wave 
function 'J (see (I0.5|) ). The corresponding ground state energy is obtained by minimizing 
this functional amongst antisymmetric TV-body wave functions : 

Ef M (TV) = inf {sf M [*], / |*| 2 = 1, * antisymmetric } . (1.20) 

I JR 3N ) 

We recall that an TV-body wave function is antisymmetric if 

i&(xi,...,Xi,...,Xj,...,x N ) = -^(xi,...,Xj,...,Xi,...,x N ) Vi ^ j. (1-21) 

For the TV-polaron there is a competition between the electronic repulsion and the effec- 
tive polaronic interaction -F C rys- There does not always exist minimizers to the variational 
problem (jl.20|) . Indeed, Frank, Lieb, Seiringer and Thomas have proved in [21 [15] that when 
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£m < 1 + a (with a independent of N), £f M has no bound state when N > 2. The following 
theorem deals with the question of existence. 

Theorem 1.2 (Ground states of Pekar's functional, N > 2, [25]). 

We assume that em > 1 in the sense of symmetric matrices. The following assertions are 
equivalent 

(1) (Binding). One has 

E^ (N) < E^ (N-k) + Ef M (k)forallk = l...N-l. (1.22) 

(2) (Convergence of minimizing sequences). All the minimizing sequences for Ef (N) 
are precompact in H 1 (M. 3N ), up to a translation. In particular, there exists a minimizer 
^em / or ^LWi solving the many-body nonlinear Schrodinger equation 

\j=i v 7 i<k<i<N 1 K J 

(1.23) 

Furthermore, for every fixed N there exists a constant cn < oo such that (|1.22j) is satisfied 
for all eu > cat. 

The above theorem is a particular case of Theorem 25 in [25] (see also Theorem 28 and 
Remark 15 therein). Frank, Lieb and Seiringer have shown recently [13] that there is also a 
minimizer if £m is the limit of dielectric matrices e„ — > £m for which til -22[) is valid for all 
n. The binding inequality (|1.22l) is a necessary condition for the compactness of minimizing 
sequences, but not for the existence of minimizers. 



1.3. Derivation of Pekar's N -polaron 

In this section we state our main results on the convergence of the polaron in the macroscopic 
limit. As motivated in the introduction, we define the energy functional for an electron 
interacting with a microscopic quantum crystal by 

£m[1>\-= \ I |VV(z)| 2 dx+m- 1 f ^ p cr (x/m)|^(x)| 2 dx+m- 1 F crys [m 3 |^(m.)| 2 ], (1.24) 

1 JR 3 JR 3 

with Vp Cr and F CTys as in Section 11.11 The corresponding ground state energy is 

E m (l) = inf |f m [^], |^| 2 = 1 j • (1-25) 
Similarly, for N > 2 the iV-polaron energy is defined as 



[*]= / Ux>^(^---> 

JR3N 



M2 , |*(a;i,...,XAr)| 2 | 

x N )\ + } : : \dxf-dx 

* — ' \Xh — Xf\ 1 

l<k<l<N 1 1 



N 



+ m 1 V°(x/m)pm{x)dx + m 1 F CTys [m 3 \p 9 (m-)\ 2 ], (1.26) 
Jr 3 

and the associated ground state energy is 

E m (N) := inf \ £ m [V], / |*| 2 = 1, * antisymmetric I . (1.27) 
I Jr 3N J 
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Similarly to Theorems 11.11 and \1.2\ it can be shown [26] that there always exists at least 
one minimizer to the small polaron variational problem ()1.25|) . and that binding inequalities 
imply the existence of a ground state for the 7V-polaron problem (|1.27[) . Here we do not need 
this and we simply introduce the concept of approximate minimizers. 

Definition 1.1 (Sequence of approximate minimizers). 

Let N > 1. We say that a sequence ($ m ) of L 2 (M. 3N ) -normalized antisymmetric N-body 
wave functions is a sequence of approximate minimizers for E m (N) if 

lim (e m [V m ]-E m (N))=0. (1.28) 
m— >0 

In order to properly state our main result, we need to introduce an Jz? -periodic function 
uf£ r which will describe the fast oscillations of our polaron, at the microscopic scale. 

Definition 1.2 (The function uf£ T ). 

We denote by uf^ r the unique positive solution of 

= inf : v E H^T), J \v\ 2 = |T| j = C'KI 



with 

1 

2m 



01 = / ^\Vv\ 2 + V p °Jv\ 2 . (1.29) 



By periodicity uf£ T is extended to the whole of M 3 . 

In Section [2] below, we will show using first-order perturbation theory that u^ r — > 1 in 
L°°(R 3 ) when m -)• 0, and that 

E^ := lim m-^r = lim m-^K"] = / (1-30) 
The function f pcl £ L°°(R 3 ) is the unique „S?-periodic solution to 

(&f pci = 2V° m . 

Ur/ por = °- 
It appears in the perturbative expansion of it*Jf r : 

||< CI '-l-m/P c l LOO(R3) <Cm 2 . (1.31) 

The following is our main result on the behavior of approximate minimizers for E m (N) in 
the macroscopic limit m — > 0. 

Theorem 1.3 (Derivation of Pekar's iV-polaron, N > 1). 

We denote by 6m > 1 the electronic dielectric matrix which was derived in flOf - Let N > 1 
be any positive integer. 

• (Energy asymptotics) We have 

(1.32) 



lim E m (N) = N £P° r + E? (AO 



m->0 



where we recall that E peT is the periodic ground state energy defined in (|1.30[) . and 



N > 2 



Ef (N) is Pekar's ground state energy defined in ()1 . 17|1 for N — 1 and m (|1.20p /or 
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• (Convergence of states) Let (\& m ) m be a sequence of approximate minimizers for 
E m {N), in the sense of Definition \l.t\ We define by the relation 



N 

y m (x u ...,x N ) = Y[u% T (x j /m) ^\x u ...,x N ). 
3=1 



(1.33) 



Then (^"'jm is a minimizing sequence for Pekar's variational problem Ef (N). 

If N = 1 or if N > 2 and the binding inequality (jl.22l) is satisfied, there exists a sequence 

of translations (r m ) C M 3 and a minimizer ^ M of Ef (N) such that 

*% i (x 1 -T m ,...,x N -T m )^*l tt (x li ...,x f ,) strongly in H\M. 3 ) N (1.34) 

along a subsequence when m — > 0. 

Remark 1.3. For this result, the fermionic nature of the electrons is not essential. The 
same theorem holds if the wave function ^ is supposed to be symmetric, i. e. if the electrons 
are replaced by bosons. 

From (|1.33j) . we see that the polaronic wave function ^ m has a multiscale behavior. The 
state is at the largest scale described by Pekar's polaron function vf^ which only depends on 
the macroscopic dielectric constant em of the crystal. It also has a fast oscillatory behavior 
encoded in the (mJz? )-periodic term uf£ T (x/m). According to (jl.31j) . this factor tends to 1 in 
L°°(1R 3 ) when m — > 0, but it contributes to V\l/ m and to the kinetic energy, yielding the term 
N E per in the total energy. The function u^ r depends on the detailed microscopic structure 
of the crystal via the electrostatic periodic potential V® eT or, equivalently, the nuclear and 
electronic densities /^p C1 . and Pp Cr . 

In the theory of homogenization, one often multiplies the microscopic periodic potential 
by m~ 2 instead of our to -1 (see, e.g., [3l|46]). The physical motivations leading to the factor 
mT 1 in front of V® cx automatically place us in a perturbative regime for the periodic equa- 
tion (|0T3)) . Actually, we could replace uf^ r (x/m) in (| X . 33|) by its first-order approximation 
1 + m f pel (x/m) without changing the result. 

For this reason, we need not use elaborate tools to analyze the contribution of the fast 
oscillations to the ground state. A simple energy decoupling method allows to separate the 
energy into the contribution of the oscillations and a functional where the fast potential 
Vp° er is absent but appears as a weight. The main point is that no derivatives of u v ^ 
appear in the latter functional. Then ()1.31|) ensures that one can ignore the weight uf£ T and 
simply consider a small polaron functional where the potential Vp° or has disappeared, i.e. the 
two scales of the problem completely decouple. The details of this procedure are provided 
in Section [2] 

The main difficulty of our work is to deal with the highly nonlinear functional -F cr y S 
which is at the origin of the dielectric matrix em and of the possible binding of polarons. 
Our main contribution in this direction is the following result, in which we show that, in a 
macroscopic regime, F crys can be replaced by under fairly general assumptions. 

Theorem 1.4 (Macroscopic behavior of F crys ). 

Let tp = (ipm) be a bounded sequence in H S {M?) for some s > 1/4. Then, we have 

Jim (TO^Fc^m 3 !^™-)! 2 ] ~ < [IV>m| 2 ]) = 0. (1.35) 
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It is rather easy to deduce our main Theorem 11.31 from this result (with s = 1), once the 
fast oscillations have been removed from the energy. The argument is provided in Section [4] 
below. 

Theorem 11.41 is based on a perturbative expansion of the energy in the macroscopic 
limit. In the case where ip m = ip is a hxed function, perturbation theory was used in |10) to 
compute the reaction of the Fermi sea of the reduced Hartree-Fock crystal on first order (this 
is the content of Theorem 3.1 of [TO], recalled above). The corresponding leading order of the 
energy was not computed however, the missing ingredient being an asymptotic expression of 
the kinetic energy that we provide in Section [3T2l below. As we will see, (jl.35|) easily follows 
from this calculation when ip m = ip. 

When ip m depends on m, then (jl.35[) is much more subtle. The Pekar interaction Ff M is an 
electrostatic interaction (similar to a simple Coulomb term) and, by the Hardy-Littlewood- 
Sobolev inequality [30 , it can be controlled by the L 12 / 5 (M 3 ) norm of ip m . By (|1.9p . the 
same is true for the crystal interaction -F C rys- It is therefore natural to assume that {ip m ) is 
bounded in L 12 / 5 (R 3 ). However, for the macroscopic limit (|1.35[) to be true, it is also crucial 
that ipm stays spread over a macroscopic region. The sequence {ip m ) could be bounded in 
L 12 / 5 (R 3 ) and still concentrate on a scale of order m. A counter example to ()1.35j) when 
(ip m ) is only bounded in L 12 / 5 (R 3 ) is provided later in Section[3?3l Eq. (|3.43|) . 

The role of our assumption that {ipm) is bounded in H S {M. 3 ), is precisely to ensure that 
{ipm) stays locally compact in L 12 / 5 (R 3 ) and does not blow up at a scale of order m, by the 
Sobolev inequality and the Rellich-Kondrachov theorem. The constraint s > 1/4 is necessary 
to have that 12/5 is subcritical, that is, below the Sobolev exponent 6/(3 — 2s). 

The proof of Theorem 1 1 .41 being the technical core of our paper, it is worth outlining its 
main steps : 

• In Section 13.21 we show by perturbation theory that one can replace the exact density 
matrix 7 m = 7p er + Q m of the crystal by its first order approximation, without changing 
-F crys to leading order. More precisely we write 

Qm — Qi + R2 

where Q\ — 0{m) and R2 — 0(m 3 / 2 ) in an appropriate norm, allowing to take into 
account only Qi when computing the leading order of the energy. This leads to an 
auxiliary functional F aux where the effect of the polarizable medium is expressed through 
a complicated but explicit operator K, studied in [lOj and defined below. This step 
relies on the application of the Cauchy formula and a resolvent expansion to (|1.2[) . 
complemented with quite a bit of algebra using the Bloch-wave decomposition of ifp Cr . 

• We then prove in Section 1331 that (|1.35[) holds, when ip m — > ip strongly in L 12 / 5 (R 3 ). 
Using the perturbative result of the previous step, this is essentially an application of 
Theorem 3.1 in [TO]. 

• We finally consider the general case. In order to quantify the intuitive fact that a se- 
quence {ipm) which is bounded in H S {R 3 ) has to stay spread over a macroscopic region, 
we use a so-called bubble decomposition. This means that we split ip m into a sum of 
pieces receding from each other in space and converging strongly in L 12 / 5 (R 3 ), plus 
a rest that is small in L 12 / 5 (R 3 ), using standard tools of nonlinear analysis as devel- 
oped by Lieb [29], Lions j33[ [34] and others. The goal is now to show that both terms 
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m~ 1 F crys [m 3 |i/) m (m-)| 2 ] and Ff u [\ip m \ 2 ] are equal to the sum of the energies of each 
piece, up to a controlled error. Since each piece of mass converges strongly in L 12 / 5 (R 3 ), 
the previous step then implies that the difference between these two sums is small when 
m — > 0. This technical argument is detailed in Sections 13.31 and 13.41 

Clearly, Theorem (ll.4[) cannot follow from simple weak limit arguments because of the 
translation-invariance properties of -Fcrys and F^ M . The term involving -F C rys does not change 
if ipm is replaced by i/) m (- — T m ) with r m in the rescaled lattice mJzf, whereas F^ M is itself 
fully translation- invariant. In particular, one can easily construct a sequence of functions 
such that VVre — "0 ( e -g- m H 1 (R 3 )), but with 

m— >0 

The natural counter-example is of the form ip m = ip + ip(. — r m ) where ip and if are fixed 
compactly supported functions and r m is a translation, |r m | — > oo when m — > 0. It is thus 
clear that our method of proof must in particular accommodate such type of behaviors and 
this is precisely what the bubble decomposition does. 

Remark 1.4. The sign of the macroscopic density v m = \ip m \ 2 > considered in The- 
orem \L^\ was only motivated by our application (the polaron). As can be seen from our 
method of proof the same result remains true if (\ip m \ 2 ) is replaced by a (real-valued) se- 
quence {v m ) with no particular sign, but with appropriate Sobolev bounds ensuring strong 
local compactness in L 6 / 5 (IR 3 ). 

The rest of the paper is devoted to the proof of our main results, Section [2] dealing 
with the periodic problem and the associated energy decoupling, Section [3] with the proof 
of Theorem 11.41 The final steps in the proof of Theorem 11.31 are presented in Section [4] 

Notation : In the whole paper, C denotes a generic positive constant whose value may 
change from line to line. 

2. The microscopic oscillations and scale decoupling 

As we have explained, our polaron has a multiscale structure. The macroscopic scale is 
described by Pekar's theory, which only depends on the macroscopic dielectric tensor £m- 
The microscopic oscillations are described by a simple linear eigenvalue problem which we 
study in this section. We also explain here how to decouple these two scales. 

We introduce an energy functional where u^ 1 (that we extend to the whole of M 3 by 
periodicity) appears as a weight. In the case N = 1, it is 

4W = \f K r (^/m)| 2 \V^(x)\ 2 dx + mr 1 F ayB [m 8 |<»| a |VK)| 2 ] . (2.1) 
In the case iV > 2 we introduce 

N 

U^(x 1 ...,x N ):=l[u^(x i ) (2.2) 

i=l 



Derivation of Pekar's Polarons 



17 



and the corresponding weighted functional is now 

(2.3) 

with the modified density 

fa(x) := K er (z/m)| 2 / TTK 1 ^)! 2 !*^^,...,^)! 2 ^...^. (2.4) 

J R 3(]V-1) ^ 

The link between these energy functionals and our original problem will be made clearer 
in Lemma \2. 2 1 

As mentioned above, the two scales appearing in our problem, namely that of the crystal 
lattice and that of the minimizer of the generalized Pekar functional decouple. Mathemat- 
ically this means that the factors involving u^ 1 in the "tilde" functional introduced above 
can be ignored because u p ^ r « 1 on the macroscopic scale. This is demonstrated in the 

Lemma 2.1 (The periodic problem). 

Let / pcr be the solution of 

W per = -V p ° cl . m r 

(2.5) 




with periodic boundary conditions on dT. There holds 

||< r -l-m/P c 1 LOO(r) <Cm 2 (2.6) 

and 

E^ = m^V° c J^ + 0(m 2 ). (2.7) 

Proof. The arguments here are standard so the the following is voluntarily kept at the level 
of a sketch. We first remark that V p ° ei € L P (T) at least for any 1 < p < 3 since it is defined 
as — AV^ cr = 47T (/ip el — Pp el ) with fi pel a non-negative finite measure, and p per £ L 2 (T). 
On the other hand, simple upper and lower bounds yield £^ er = 0{m). From the equation 

--Au pcr + mV° u pcl = mE pel u pel 

we then deduce that ||uP^ r — l|| ff i(n = 0(m). Writing u p ° r = 1 + mf it is then easy to see 
that necessarily f r f = 0(m) and 

-±Af = E p ?-V p CI + O LP (m) 

for any 1 < p < 2. This implies that ||/ - / pcr || L oo (r) < Cm (recall that E^ 1 = 0(m)) by 
standard elliptic regularity and the conclusion of the Lemma follows. □ 



The following lemma is the key tool allowing us to effectively decouple the two different 
scales of the problem. Such decoupling techniques have already been used in several contexts, 
e.g. Ginzburg-Landau theory |23j and the derivation of Gross-Pitaevskii theory from many- 
body quantum physics [31] , 
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Lemma 2.2 (Energy Decoupling). 

Let ^ 6 H 1 ^ 31 ^) be a normalized wave function, J m3N \^\ 2 = 1- Define V I' P 6y i/ie formula 

* = E/r(-/m) * P (2-8) 
where uf" is extended to R 3 by periodicity and Uf° r {x\, xjy) = rij=i ^m T { x i)- We have 

S m [*] = Nm^E^ + S^}. (2.9) 

Proof. Let us first remark that ^ p is well defined by formula (|2.8p , because is strictly 
positive (it is the ground state of a Schrodinger operator). 

We only write the proof for N = 1. Using that = Pu^{-lm)^i the argument is exactly 
the same for N > 2. The only input of the lemma is the variational equation satisfied by 

- iA< r + m^C = m^ST. (2-10) 
From the definition (12.81) and an integration by parts, we have 

f |V*| 2 = m- 2 / (-A< or (./m))< or (./m)|* p | 2 + / !<"(•/»») | a | W p | 2 . 

JR3 JR 3 JR 3 

Hence, using (|2.10[) and the normalization of "J, we obtain 
\f iW&lVm- 1 / U p ° cr (7m)|*| 2 +F crys [m 3 |vI/(m.)| 2 ] 

Z JR 3 JR 3 

= m- i £;r+ / |<r(-M)| 2 |V* P | 2 +F crys [m 3 |un 2 |* P K)| 2 ]- 

,/R 3 

This gives the result in the case N = 1. □ 
3. Derivation of Pekar's interaction: proof of Theorem 11.41 

In this section we provide the proof of Theorem ll.4[ which is our main result relating, in the 

macroscopic limit, the exact crystal interaction i^crys to the generalized Pekar interaction 
F p 

em 

We start by recalling in Section 13.11 several properties of the crystal model that we use 
in our analysis. In particular we introduce an auxiliary model that is 'halfway' between the 
full interaction via the crystal and the Pekar macroscopic interaction. In Section 13.21 we 
explain how to go from the crystal model to the auxiliary model. The latter is a Coulomb- 
type interaction involving a complicated but explicit operator denoted K,. This operator K, 
is nothing else but the linear response of the crystal in the sense that, if pq denotes the 
perturbation of the crystal density in presence of the defect v, one has pq ~ — TCv for v 
small enough (see (|3.18|) below). Results about JC proved in [10] (alluded to briefly in the 
introduction and recalled below) suggest that the intermediary model naturally turns into 
a Pekar theory in the limit. Namely, they allow, in combination with the analysis of Section 
13.21 to prove (|1.12|) . However, the proof of Theorem 11.41 requires to implement these ideas 
under much more general assumptions than was done in |10j . This is done by using a bubble 
decomposition for the sequence (4>m), which is the content of Sections 13.31 and 13.41 
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3.1. Useful properties of the crystal 

In order to realize our program, we first recall several useful facts about the polarized crystal 
in an external density v, whose proofs can be found in [8| 110] . 



3.1.1. Linear response of the Fermi sea 
We denote by 

^"crysIQ, V] = D(l>, p Q ) + T CIys [Q] (3.1) 

the total energy of the crystal in presence of the defect v, with J- C ry S [Q] as in (11.7[) . and we 
recall that 

■FcrysM = hlf ^"cry S [Q,H- ( 3 - 2 ) 

-7p Cr <Q<l-7por 



The proper functional setting for the minimization problem Q3.2p is detailed in [8], we only 
sketch it here. Any operator Q satisfying the constraint 

" 7pcr < Q < 1 - 7pcr (3-3) 

is decomposed as 

Q = Q— +Q- + + Q ++ + Q + - (3.4) 

where Q~~ = 7p Cr ( 97pcr> Q h = IpaQ (l ~ 7pcr)j an d so on. It is proved in [8] that for 
Q satisfying ([331) and v <E L^R 3 ) n L 2 (R 3 ), T ay& \Q, v] is finite if and only if Q is in the 
function space 

Q = [Q g & 2 \Q = Q*, |V|Q G 6 2 , Q ++ ,Q— G 6 1 , |V|Q ++ |V|, |V|Q— |V| G S 1 } . 

(3.5) 

The space Q is endowed with its natural norm 

IIQIIc = HOIIe- + IIQ ++ ll el + HQ""llei + ll|V|Q|| e 2 + |||V|Q++|V||| el + |||V|Q— |V||| el . 

The symbols 6 1 and & 2 denote the Schatten classes of trace-class and Hilbert-Schmidt 
operators on L 2 (R 3 ) respectively (see [H] and 03]). 

An important fact is that, although the operators in Q are not necessarily trace-class, 
they are always locally trace-class and thus have an unambiguously defined associated den- 
sity pq G Lj^R 3 ). In Proposition 1 of [8], it is even shown that pq belongs to L 2 (R 3 ) and 
to the Coulomb space 

C= {p\D(p,p)^ 2 <co}. (3.6) 

By definition, there holds 

Tro(^Q) - f V PQ (3.7) 
for any V G C (the generalized trace Tro is defined in (|1.6p ). Moreover the inequality 

WpqWl^+D^Pq) 1 ' 2 <C\\Q\\q (3-8) 



holds uniformly for Q G Q, showing that the linear map Q G Q n- pq G L (R ) n C is 
continuous. 
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The existence of Q £ Q minimizing Q3.2p for fixed v is proved in [8], Theorem 2. For ease 
of notation we will not emphasize the dependence of the minimizer of (|3.2[) on v. Minimizers 
are found to satisfy the self-consistent equation 

Q = 1 ( _ 00>£F) (ff° er + {pq + v)*\. - Ic-oo^) (# p ° or ) + ^ (3.9) 

where 8 is a finite rank self-adjoint operator with 

Ran(5) c Ker (F° er + (p Q + v) * \ . I" 1 ) . 

On the other hand, there holds 

> -F CI ys[Q, V\ > \d{ Pq , Pq ) + D{V, p Q ) > D( PQ) pq) 1 / 2 ( \d( P q, pq) 1 ' 2 - D(v, v) 1 ' 2 



where the upper bound is obtained by using a trial state Q = and the lower bound by 
dropping the positive kinetic energy term and using the Cauchy-Schwarz inequality We 
deduce the uniform estimate 

D{p Q ,p Q )<AD{v,v). (3.10) 

An easy consequence of the bound (|3.10[) is that the difference between H® cr . and H® er + 
(pq + v) * | . | _1 is controlled by D(v, v), in the resolvent sense. In particular, if D(y, v) is 
small enough (depending on ep, see [10] Lemma 5 for the short proof) then Tro(Q) = and 
e F i o-(ff° er + (pQ + v)*\. I' 1 ). Then 8 = in ([379]) and Q is unique: 

Q = l(-oo, £P ) (^per + (PQ + f) * I ■ I"') ~ l(-ao,e P ) (^pcr) ■ (3-H) 

Since £f is not in the spectrum of the mean-field operator when D(i>, v) is small enough, 
it is possible to express Q using Cauchy's formula. Let ^ be some curve in the complex 
plane that encloses the spectra of both i?p er and H® er + (pq + v) * | . | _1 in (-co, Using 
the resolvent formula one can then decompose Q as 

Q = Q 1 + R 2 (3.12) 

with 

Qi:= i£^^ y ^W R2 := i£z-i-v( v i^) (3 - i3) 

and 

V :=(p Q + J /)*|.|-i. (3.14) 

We again refer to [TO], in particular Section 6.3, for a more precise discussion of these facts. 
When D(y,v) is small, V is also small in an appropriate sense because of f|3 . 10|) . One can 
then prove (see the proof of Lemma 3 in [TO]) that 

IIQi|| c <c||y|| L2+c , <CD{ v , v yi 2 , 

\\R2\\ Q <C\\V\\ 2 L2+c ,<CD(v lV ), (3.15) 

for all D(v, v) small enough. 

The main idea is now to just discard the second-order term R 2 and only keep Qi in our 
energy functional, leading to the auxiliary interaction F^^. This is done as follows. First we 
consider the first order linear operator 

1/1 1 

:= -PQiW, where Qi{p) := — — (p*\ ■ \~ ) — ■ (3.16) 
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It was proved in |10j (Proposition 2) that £ is a bounded non-negative self-adjoint operator 
on the Hilbert space C . Explicit formulas can be found for the operator £ in [10] but they are 
not needed here. Assuming that D(is, v) is small enough, we can rewrite the self-consistent 
equation as follows 



PQ 



-£(pq 



PR2 



or, equivalently, 



PQ = -(l-{l + C)- 1 )v+{l+C)- 1 



PR2 



Introducing the operator 



/C:=l-(l + £)- 1 . 



(3.17) 



which is bounded and non- negative on C, we deduce from (|3.8I) and (|3.15[) that 

\p Q +Kv\ c < \\p R2 \\ c <C\\R 2 \\ Q < CB{v,v) (3.18) 

for D(y,v) small enough. The leading order for pq is therefore —JCv, which is interpreted 
as the linear response of the nonlinear Fermi sea. 

We now define the auxiliary interaction -F a ux by retaining only second order terms in the 
energy, which corresponds to keeping only the first order terms for Q (the first order energy 
automatically vanishes). The correct formula turns out to be 



■-D(v,Kv), 



(3.19) 



as we will prove in Section [3.21 below. 



3.1.2. Macroscopic dielectric matrix 

Other results of [TU] will prove useful in our context. We define the dilation operator U m 
and its adjoint in L 2 (R 3 ) by 

{U m v){x) = m 3 v(mx), (U^v){x) = u(x/m). (3.20) 

We then introduce the potential 

B m {v) = m- x U: n (| . T 1 * (1 + Cy 1 {U m i/)) . (3.21) 

The above formula means that given a charge-defect i>, we first scale it using U m , apply the 
operator (1 + C) , consider the corresponding electrostatic potential (convolution with the 
Coulomb kernel) and then scale back the total result using U^. An easy computation shows 
that 



m 1 F aux [m 3 i/(m-)] 



B m (v)v - D(u, v) 



(3.22) 



Using that C > as an operator in the Hilbert space C (Proposition 2 in [10]). implying 
that (1 + £) _1 < 1, we obtain 

0< J B m {v)v = m- 1 D(U m v,(l+L)- l U m v)<m- 1 D{U m vM m v)=D{v,v). (3.23) 
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This proves that B m is continuous from C to C, with a norm independent of m: 

\Bm\ c ^c> < !• (3-24) 

Moreover, from the analysis in Section 6.10 of [10], we have for any fixed v € C that 

B m (y) -» W v (3.25) 

weakly in C, where W v is the dielectric potential depending on the dielectric matrix em, 
which we have defined above in (|1.15[) by 

- V (emVW^) = Attis. (3.26) 

This result is a step of the proof of Theorem 3.1 in [10], whose main conclusion was recalled 
in Section ITTTl 

The connection between the full interaction with the crystal and Pekar's model goes 
roughly speaking as follows : 

mT 1 T crys [m 3 'v{m.)\ s» m _1 F aux [m 3 i/(m.)] 
1 



B m {v)v-D{v, v) 
\(JjV v v-D{v,v)}=Fljy\. 



Section 13.21 is devoted to the proof the first 'equality'. The second equality is a simple 
scaling. The convergence (13.25)) then proves that the third 'equality' is justified when v does 
not depend on m. To allow the m dependence specified in the assumptions of Theorem ll.4l 
we need other ingredients, that are provided in Section 13751 and [3T4l 



3.2. Perturbation theory and the auxiliary interaction F aux 

In this section we use perturbation theory to approximate the nonlinear reaction of the crys- 
tal (described by the operator Q), by its first order approximation, obtained by considering 
that the defect v is small. This will allow us to replace the complicated nonlinear interaction 
energy F QTys by its leading order -Fa UX . For completeness, we consider an abstract defect v 
(for instance in L 6 / 5 (R 3 )) and we provide uniform error bounds. 

We will come back later to the case of a sequence v m = m 3 \ip m (m-)\ 2 with (ip m ) bounded 
in iP(R 3 ), as is considered in the statement of Theorem ll.4l As was already explained and 
used in [10j . the correct norm measuring the size of a given defect v is that given by the 
Coulomb term D{v, v). Saying differently, one can prove that the terms in the perturbation 
series of pq can all be controlled by powers of D(i>, v). Fortunately, we have 

)| 2 ,m 3 |^ m (m-)| 2 ) =mD{\^ m \ 2 ,\^ m \ 2 ) < Cm ||^ m ||£ 12/5(1{3) < Cm, (3.27) 

when (ip m ) is bounded in H S (U. 3 ) with s > 1/4. The macroscopic limit m — > thus auto- 
matically places us in this perturbative regime, which is why we will be able to only keep 
the leading order term. 

The following proposition justifies the introduction of the auxiliary interaction functional 
by quantifying its difference with the full interaction with the crystal. 

Proposition 3.1 (From F crys to F aux via perturbation theory). 

There exist two positive constants C and rj such that for any v G C with D{y^v) < r\, we 
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have 



-^crys [ 



<CD{v, vf' 2 



where F &ux was defined above in (13.191) . 



(3.28) 



The rest of this section is devoted to the proof of Proposition 13.11 In the first step we 
rewrite the kinetic energy following [5J [2U] . 

Lemma 3.1 (Rewriting the kinetic energy term). 

Under the assumptions of Proposition [PI we have 

Q ++ - Q~ = Q 2 , (3.29) 

and consequently 

Tr ((ff p ° cr - ep) Q) = Tr (\H° pcl . - s F \ Q 2 ) . (3.30) 
Proof. Recall that there holds 

Q = l ( _ oc , £F) {H° pcI + (pQ+l>)*\. I -1 ) - t(-oo,e F ) K°er) 

when D(i/, v) is small enough. We write this as Q = II2 — III for short, with obvious 
definitions. Then Q 2 = II2 + III — IIin 2 — n 2 IIi using U 2 = Hi. But in this notation we 
have Q— = nin 2 ni - III, and Q++ = U 2 + nin 2 ni - n 2 rii - nin 2 which yield (|339]) . 
Then (|3.30[) is obvious from the definition of the kinetic energy term Tro ((-ffp Cr — £f) Q) — 
Tr(|i/ p ° or - £F |(Q++-Q-)). ' □ 

The second step is less easy. It uses the decomposition Q = Q\ +i? 2 and the special form 
of Q\. As we explain below, the result should be true in a more general framework, but we 
present its proof in our context only, using the Bloch decomposition of H pcr and the fact 
that the spectrum of this operator has a gap as crucial ingredients. 

Lemma 3.2 (Leading order of the kinetic energy). 

There holds 

Tr (\H° pm . - e F | Q 2 ) = ~ Tr (QiV) (3.31) 
where V and Qi are defined in (|3.14[) . 

To make the result plausible, let us present first a calculation which is valid for matrices. 
Let C be an hermitian matrix and let us denote 

nc = i ( -oo,o)(C)- 

Writing for any operator Q 

Tr(CQ) = Tr (|C| ((1 - II C ) Q (1 - II C ) - II c Qn c )) 

it follows that 

Tr (C(H - n c )) = Tr (|C|(II - n c ) 2 ) , (3.32) 

for any operator C and projector II. The proof goes along the very same lines as that of 
Lemma [33] with the fact that II = II 2 as only input. 



2-1 



Let us now consider two operators A, B and denote 

Q = ^A+mB — ^A- 

We will assume that m € R is a small parameter and that Q can be expanded as 

Q = mQ + 0{m 2 ). 

In the context of Lemma 13. 2\ 

A = H° CI -£ F , m = D( Vl v) 1/2 , B = m~ 1 V, Q = Q x /m 

with V defined in (|3~T4| and Qi in ([37T3]) . Using ([332]) with C = A + mB and n = U A we 
first compute 

Tr {(A + mS)Q) = - Tr (\A + mB\Q 2 ) . 
Inserting the expansion of Q we thus have (at least formally) 

Tr ((A + mB)Q) = - Tr (\A + mB\Q 2 ) = -m 2 Tr(|,4|Q 2 ) + (9(m 3 ). (3.33) 
On the other hand, by linearity and using (I3.32[) with C = A and II = TLA+mB we have 
Tr ((A + mB)Q) = Tr (AQ) + m Tr (BQ) 

= Tr (|A|(3 2 ) + m Tr (BQ) 

= m 2 Tr (|^1|Q 2 ) + m 2 Tr (bQ \ + 0(m 3 ). (3.34) 
Comparing (|3.33|) and (|3.34j) we obtain, to leading order in m 

Tr(\A\Q 2 ) =~Tr (BQ) 

which is (j3~3T|) . 

It is elementary to justify the previous calculations for finite matrices, when is not in 
the spectrum of A, that is when A has a gap in its spectrum |22j . In the infinite dimensional 
setting of this paper, the operators are not bounded and the trace Tr has to be replaced 
by the generalized trace Tro- This generates serious difficulties if one tried to put the above 
sketch on a rigorous basis. For this reason we do not follow the simple strategy presented 
above and provide instead a proof of Lemma 13.21 relying on computations in a Bloch basis 
diagonalizing A — H® CI — £p-- 

Proof. We use the Bloch- wave decomposition of -ffp Cr [43] : 

=-f ((H° pcl ) q f q ) e^dq. (3.35) 

with 

(H° eI ) q (x) = J2 V? K ? > (a n J ■ (3.36) 

Then we have (recall (jl.ip ) 

+oo 

(Tper)?^) = Yl 1 (-oo, £P )(An ) g) K, ? ) (^,9! 
n=l 
+ 00 

(! - Ipcr)q(x) = l( £F ,+oc)(A„, g ) |a„, ? ) (a„, g | . (3.37) 

n=l 
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We remark that as V is a multiplication by a potential there holds [V] p , q = (V) p - q which 
allows to compute the components of the Bloch matrix of Qi, starting from (|3.13[) 

= e £ (,-aJz-vJ (a <-' |a - } (a -' 

= X] 1 (eF,+oo)(Ai,p)l(-oo, £F )(A„, 9 )- — (o/, p | (V)p-g |On,g) |o/,p) (ffln.gl 

l.n A ™'« 

+ X 1 (-° - e F)(Az,p)l(£ F , + oo)(A„, ( j)- ZT\ ^'pl l a ™^) Kf) ( a ™'9l 

= [CU + [D] p , q . (3.39) 

We have used the residuum formula to evaluate the integral over the curve . We then note 
that the operators C and D defined in the above formula satisfy Ran(C) C Ker(C) and 
Ran(D) C Ker(D), so Q\ = CD + DC. 

We compute the components of the Bloch matrix of CD using the formula (the symbol 
f = irT 1 / r . denotes the mean value over the Brillouin zone T*) 

[CD]p^ q — + Cp^pi Dp/ 19 
Jp' 

and obtain 

[CD] Pj g = 2J l(£ F ,+00)(A/ : p)l(-00,e F )(A n ,p') 



, ti , (A n , P ' Aj )P )(A n ,p' A n ' j9 ) 

x (ai, p \ (V)p-p' |a«,p') (a„,p'| (V)p/-? | 
Now [|i?° er - e F | CD] = |iJ° er - £ F | p [C-D]p, P and thus, recalling ([335]) and (j3~36l) . 



Tr (\H% -e F \CD)=-f Tr (|ff° er - e F | [CD] p , 3 



X 1 (£ F ,+oo)(A;,p)l(- 00 , £F )(An,p') — — -^—2 
P' i, n (A„ :P ' - Ai,p) 

2 



x l( a *,pl (^)p-p' l«n,p')T 



where we have used that (a;^] (V)p_p/ |a n ,p') = ( a n,p'| (^)p'-p l a z.p)- A similar computation 
leads to 

Tr (l-ffper - £ f| DC) = -f -f l( EF ,+ 00 )(A;,p)l(- 00 , £F )(A n ,p') " ,y — ^-5 

P P' in (.Aj,p — ^n,p'J 

x l( a i,pl (^)p-p' l a «,p')| 2 ■ 
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Summing the two contributions and noting that for A;. p > ep and \ n ,p' < £f we have 
|A/, P — £f| + \Ki,p' — £f| = A; lP — \n,p', we conclude that 

Tr (l^pcr ~ £ f| Q 2 ) = /t ^ 1 (e F ,+oo)(Ai 1 p)l(- o 1 e P )(A n ,p')-r 

Jp-'P' I n ' P ■ P ' 

x \{a hp \(V) p - p < |a„, p ,)| 2 . (3.40) 

An independent computation starting from (|3.38p yields 
Tro(FQ) = / / E <°''.pI (^W(Qi) 

= // ^ 1 (ep,+oo)(Aj,p')l(-oo,ep)(An,p)T — « l( a i,p'l (K)p'-p |<Vp)| 2 



•^P i,n 



l(e F ,+oo)(A„,p)l(_oo : e F )(A(,p')T — T K^pl 00p-p' l<Vp')l 

= ~ 2 / / 1 (ef,+oo)(A/ : p)1(- o, £f )(A„^)- — — \{ai, P \ (V)p-p' \a n ,p>)\ 2 

and proves the lemma. □ 

We are now ready for the 

Proof of Proposition 13.11 With Lemmas 13.11 and 13.21 at hand it essentially remains to 
bound the contribution of R2 to the kinetic energy. This is done by using the estimates 

Tr (\H% eI -e P \AB) < C\\A\\ Q \\B\\ Q . 

and 

\TMUA)\<C\\A\\ Q \\U\\ L2+C , 

that hold true for any operators A,B and any potential U(x), see [8J. Using these in com- 
bination with (|3.10p and p,15[) leads to 

|Tr (\H° eI - £f|Q 2 ) - Tr (|if p ° er - e ¥ \Q\) \ < CD{v, vf' 2 
\Tr (VQ) - Ti (VQ!)\ < CD{v, vf' 2 . 

Using that by definition of V, Tro(VQ) = D(pq 1 pq + v), and combining Lemmas 13.11 and 
13.21 we infer 

F ciys [v] = J" crys [Q, v\ = -^D(p QlPQ +v) + D{v, PQ ) + \d{ Pq , Pq ) + O (d(u, vf' 2 ) 

= l -D{p Q , V )+o(D(v,vY' 2 ). (3.41) 

Now, because of (|3.18[) . we can replace pq by — JCv in (|3.4ip at the expense of an other error 
controlled by D(v, v) 3 ! 2 . We thus obtain 

F ciys [u] = ~D(u,Ku)+o(d{u,u) 3 / 2 
as was claimed. □ 
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3.3. Towards the effective interaction F^: two lemmas 

Coming back to our main question, we have by Proposition 13.11 and the estimate (I3.27P 

m _1 .Fc^[m 3 |^ m (m-)| 2 ] = m-^^m^mK)! 2 ] + 0(m l l 2 ), (3.42) 

for any bounded sequence (ip m ) m H S (M. 3 ) with s > 1/4. It therefore remains to study the 
limit of the right side of (ET32|) . 

For a fixed defect v £ C, we deduce from p.42j) and (|3.25p that 



lim m 1 F c ^y S [rn 3 v(rn-)\ — lim — [ f 



B m {v)v - D(v, v) 



We are now interested in proving the same kind of result for an m-dependent sequence of 
the form v m — li/Vrel 2 , which is of course much more complicated. This is indeed not true in 
general, without appropriate assumptions on (ip m )- 

For instance, only assuming that (v m ) is bounded in C is certainly not sufficient to get 
the limit. Consider the example 

v m = m 1/2 m^v(x/m) (3.43) 

for a fixed density v 1 where the scaling is chosen such as to have D(v m , v m ) = D{v, v). In 
this case we have by (|3.23p 

B m {v m )v m = m~ 1 D(U m is m , (1 + Cy^mVm) = D{v, (1 + £)~M 

which is not close to 



/ W Vm v m = I 



In this example the sequence v m actually lives at the microscopic scale because of our chosen 
scaling and a macroscopic behavior cannot be expected. 

In order to be able to approximate m~ l F crys [m i v rn (m-)] by its macroscopic counterpart 
to -1 ^ ( [TO 3 ^ m (r7i-)], we have to make sure that v m stays spread over a macroscopic region 
when m — > 0. One way to ensure this is to impose that v m is locally compact in L 6 / 5 (IR 3 ). 
For v m — IV'ml 2 , this follows when ip m is bounded in H S (M. 3 ) for some s > 1/4 : 

Proposition 3.2 (From F aux to Ff in the locally compact case). 

Let (ip m ) be a bounded sequence in H S (M. 3 ), s > 1/4. We have 

lim (m- 1 ^™ 3 ^™-)! 2 ] - F £ P M [|^ m | 2 ]) = 0. (3.44) 

Combining (|3.42p and (|3.44p concludes the proof of Theorem 11.41 In the context of our 
polaron model, natural a priori H 1 bounds will be satisfied by minimizing sequences, so the 
result will be applied with s = 1. 



The proof of Proposition 13.21 uses the local compactness by resorting to standard tech- 
niques of nonlinear analysis. We write ip m as a sum of pieces converging strongly (up to trans- 
lation) in L 12 / 5 (]R 3 ), and receding from each other, plus a rest that is small in L 12 / 5 (R 3 ) (a 
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method that is usually called a bubble decomposition in the literature). In order to conclude, 
we then need to prove two things : 

(1) that one can pass to the limit for each compact piece of mass 

(2) that the interaction energy between two such pieces is negligible. 

In this section we prove two lemmas going in this direction, deferring the bubble analysis in 
itself and the conclusion of the proof of Proposition 13.21 to the next section. 

Our first lemma deals with the case when v m converges strongly to some v in C, up to a 
translation. 

Lemma 3.3 (Self-interaction for strongly convergent sequences). 

Let (v m ) be a sequence of functions such that 

v m —> v strongly in C (3.45) 
when m — > 0. Let (x m ) C R 3 be any sequence of translations. There holds 

lim / B m (v m (- - x m )) v m {- - x rn ) = / W v v. (3.46) 

m-H) J R3 J R 3 

Proof. Using the fact that B m : C — > C is bounded uniformly in m by (|3.24|) . we deduce 
from the strong convergence of v m 

B m {v m {- - x m ))v m {- - x m ) = / B m (v(- - x m ))v{- - x m ) +o(l). 

The lattice spacing of the crystal is of order m in our scaling, thus it is always possible to 
write 

where r m £ m«5f is a translation of the scaled lattice and \y m \ < Cm. It is clear from 
the definition of C in [8] that C commutes with the translations of the original (unsealed) 
lattice, and this implies by our definition (|3.2ip that B m commutes with the translations of 
the scaled lattice ?tiJz? . Thus 



B m {v{- - x m ))v(- - x m ) = / B m (v(- - y m ))v{- - y m ). 
Note that — y m ) — > v strongly in C, hence we arrive at 



B m {v m (- - x m ))v m (- - x m ) = / B m (v)v + o(l) 

and there only remains to use the result (|3.25l) proved in [10] to conclude the proof. □ 

We now provide the basic result allowing to neglect the interaction of two pieces of mass 
receding from each other. Note that the Pekar interaction of two such pieces of mass is easy 
to estimate because of the explicit form of its kernel, but the auxiliary interaction F &ux has 
no such simple form and the argument is more involved. We actually have to come back to 
the technique of [10] used in the proof of (|3.25p . 

Lemma 3.4 (Interaction of two pieces of mass receding from each other). 

Let {v m ,i) and (v m ,2) be two sequences of functions such that 

v m ,i — > v% strongly in C, (3.47) 
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for i = 1, 2 w/ien m — > 0. Let (a; m) i) and (x mi 2) 6e iwo sequences of translations in R 3 such 
that \x m> \ — x m _2\ — > +oo. There holds 



lim / B m (v m l (- - x m ,i))v m 2{- - x m 2 ) = 0. (3.48) 
m->-0 ,y R 3 

Proof. As in the proof of the previous lemma, the uniform boundcdncss of B m and the 
strong convergence p.47|) yield 



B m (i* m .i(- - x m ^))v ma {- - x m ,2) = / B m {vi(- - x mt i))v2{. - x m ,i) + o(l). 

For the rest of the proof, we re- use some arguments of Lemma 7 in [10]. We define for i = 1, 2 

0»,m = w c 1/2 ^(- - z m ,i) £ L 2 (M 3 ) 
where v c is the Coulomb operator: 

v c (v) = v*\ . I" 1 . 

By density of regular functions in Fourier space, we can assume that both v\ and v-i are in 
C c yo (M 3 \ {0}). Following [10] Lemma 7 we then have, for m small enough, 



5 m (i>l(- - x mi i))v2{. - x m> 2) = I dk((e 1 ) mfc e ,e )5T^ l (A:)5^n(fc) 

dkdr 1 )^ e , e ) Wk)92{k)e^-^-\ 

with e = v c (1 + £) w^ 1 , e = i^^e Vc^ 2 and = v X J 2 Vi. In the second line we have used 
that v c commutes with multiplications in the Fourier space (it is itself such a multiplication). 
The factor ( (£ _1 ) mfc e o, eo) is uniformly bounded and tends to \k\ 2 /k T eMk a.e., as proved 
in |10j . Lemma 6. By the dominated convergence theorem and Riemann-Lebesgue lemma, 
we thus have 



lim / B m {v 1 {--x m , 1 ))v 2 {.-x m>2 )= lim / rffc ^ff } e^"^-')'* = (3.49) 
because e 1 ^" 1 ' 2-21 '"' 1 )'' converges to weakly-* in L°°. This concludes the proof. □ 



3.4. Completion of the proof of Theorem \l-4\ 

The proof of Theorem 11.41 uses some classical tools of nonlinear analysis [29] |33l [34] [35] , 
which are for instance recalled in the appendix of [24] whose notation will be used here. The 
starting point is the 'highest local mass' of a sequence u = (u m ), following ideas of Lieb [29] . 

Definition 3.1 (Highest local mass of a sequence). 

Let u = (u n ) be a bounded sequence in iJ s (K 3 ). We define 

M(u) = sup { / M 2 , 3(x k ) C U 3 : u nk (• - x k ) u 

weakly in H S {S?) for some subsequence (n^) >■ (3.50) 
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As first noted by Lions [33l [34] , the above quantity plays a central role when controlling 
subcritical LP norms of the sequence, see Lemma 8 in [24] , 

Lemma 3.5 (Control of subcritical LP norms). 

There exists a constant C , depending only on s, such that 

f 4s 2s 

limsup / |w„| 2+ ~ < CM(u)Tlimsup||u„||^ s(K3) . (3.51) 

n— foe JIR 3 n—>oo 

We now recall a classical result allowing to locate the "compact pieces of mass receding 
from each other" we were alluding to before. Such bubble decompositions are ubiquitous in 
the literature (see, e.g., [UJ [35J [IB])- The following statement is taken from fM\ . 

Lemma 3.6 (Bubble decomposition of a bounded sequence in LL S (M. 3 )). 
Let u = (u n ) n be a bounded sequence in H S (M. 3 ) with 1/4 < s < 3/2. For any e > and any 
sequence < R n — > +oo, there exists JeN, sequences of functions uj = (uj, n ) n , j = 1 • • - J 
and Uj = {Uj_ n ) n , sequences of translations Xj = (xa n ) n , j — 1 ... J such that, along some 
subsequence (we omit to change indices) 

.7 

u n — Uj >n (- — Xj >n ) — Uj.n —> strongly in LL S (M. 3 ). (3.52) 
i=i 

Moreover we have 

(1) Uj tn -> Uj - ^ weakly in LL S (R 3 ) and strongly in L P (R 3 ) for any 2 < p < 6/(3 - 2s) 

(2) supp(uj. n ) C B(0, R n ) for any j — 1 . . . J and any n e N 

(3) \xj in — Xk, n \ > 5i? n for any j ^ k and any neff 

(4) supp([/j, n ') c R 3 \ Uj =1 B(x jtn , 2R n ) 

(5) M(Uj) < e. 

We now complete the proof of Theorem 11.41 The only missing ingredient is the 
Proof of Proposition 13.21 We appljQ Lemma 13.61 to the sequence (ipm), and denote 
v m = \ipm\ 2 - We thus have, along some subsequence 

J 

\\lpm - Vm||,ff»(R3) °) witn <Prn = ^2<*P],m(- ~ Xj,m) + *J,m- (3.53) 

3=1 

The functions <Pj,m, j = 1 . . . J, VP j jm and the translations Xj >m satisfy Items (1) to (5) in 
Lemma 13.61 for some e, R m . We denote tpj t oc the limit of the sequence (pj t m- We further note 
that, as the functions appearing in (|3.53p all have disjoint supports, there holds 

,7 

\(Pm\ 2 = ^ l^ m (' ~ X 3\™)\ 2 + l*.7,™| 2 ■ (3.54) 

Notice that the assumption s > 1/4 implies that 12/5 < 6/(3 — 2s) hence we have |<£j,m| 2 — > 
I fj,oo 1 2 strongly in C, by the Hardy-Littlewood-Sobolev inequality [30] ■ Also, using (|3.51[) 
and an interpolation argument, we can control the L 12 / 5 norm of the rest ^ j, m by 

limsup|*j, TO | Lla/8(R3) < CM(*j)3+2 S limsup ||*j, m || ff s( 3 + 2s < (3.55) 



4 Strictly speaking wc first extract a subsequence so as to consider a discrete set of m's. 
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where 9 > is such that 5/12 = 30/(6 + 4s) + (1 - d)/q with q = 2 for 3/10 < s < 3/2 and 
g = 6/(3 — 2s) for 1/4 < s < 3/10. The number e is that appearing in Item (5) of Lemma 



We now claim that (along the appropriate subsequence) 



B m (|Vm| 2 )|^m| 2 -^ / W\ VjoB{ *\<p jlC 
j=l "'K 3 



(3.56) 



< Ce3+2s +o(l) TO _x). 

To this end we use the decomposition (|3 . 53|) and Q3.54p to obtain 

B m (\lpm\ 2 )\lpm\ 2 = V] f B m (\tpj, m (- ~ Xj, m )\ 2 ) \<pj,m(- ~ X^ m )\ 2 + [ B m (|*j, m | 2 ) \Vj, m \ 
~[ JR 3 JR 3 

+ /2 B m (\cpj. m (- - Xj. m )\ 2 ) \(p k , m (- - X k ,m)\ 2 

V / B m (\(p jim {- - x Jtm )\ 2 ) |*j, m | 2 + o(l). 

• -, JR. 3 



(3.57) 



Now, using the continuity of B m from C to C and the Hardy-Littlewood-Sobolev inequality, 
we obtain from (|3.55[) 

40 s 



/ B m (|*j, m | 2 ) |*j, m | 2 < OD(|*j, m | 2 , |* j, m | 2 ) < C 

JR 3 

Similarly, 



m|| i l2/5( R 3) < Ce3+2s ( 



V / B m (\ipj, m (- - x^ m )\ 2 ) \^j, m ? 



< C 



J'=l 



26s 



||*j,m|l2i a /»(R3) <C£ 3 + 2s - (3.58) 



L6/5( K 3) 

In the second line we have used f|3 . 54[) to infer that 

^ ] Wj.mi' — x j.m)\ 



- II l^ m | 2 |L<!/5(R3) < C) 



L6/5(R3) 



independently of J. On the other hand 

lim / B m (\<Pj, m {- ~ Xj,m)\ 2 ) \<Pj,m(- - %j,m)\ 2 : 
m ^°jR 3 

and, for j ^ k, 

lim / B m (\<pj, m (- - ^j. m )| 2 ) \<Pk,m(- ~ a^,m)| 2 







by a direct application of Lemmas 13.31 and 13. 4[ respectively. We have thus proved the claim 
(13361) . 

A similar (but simpler) argument shows also that 



JR 3 j=l 3 



29s 



< Ce3+2s +o(l) r 



(3.59) 
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Indeed, one can decompose L 3 Wi^, ia \ip m \ 2 like in Q3.57p . and for any v, [i we have 



-v(k)p,(k) 



< C 



j-nal^O)! IA0)I < C Ikllie/s^a) lMli6/s ( 



4ir J R 3 fc T £Mfc 

All the terms can then be dealt with exactly as before. 

As a conclusion, comparing (|3.56l) and (I3.59|) . we have shown that 



lim sup 

m— iO 



B m (\^m\ 2 M 



26s 



Since e can be chosen as small as we want, this proves the statement (for a convenient 
subsequence, but standard arguments yield the result along the whole sequence). □ 



4. Proof of Theorem [T31 

We have now all the tools needed for proving our main result. Let (^ m ) be a sequence of 
approximate minimizers in the sense of Definition 11.11 We write 



* m = Ur(-/m)^ 
and use the energy decoupling of Lemma 12.21 to obtain 

£m[*m] - Nm- l E%* +4K']- 
As a consequence of the L°° estimate (|2.6[) we deduce 

£ m [* m ] = Nm- l E%* + ^[^(l + 0{m)) 

with 

l*l 2 



(4.1) 



(4.2) 



E 



i<k<e<N 



3N \Xk - Xt\ 



+ m 1 F crys [m 3 pq,(m-)] 



(4.3) 

The meaning of ()4.2p is that the polaronic problem encoded in the functional (14. 3|) has been 
completely decoupled from the microscopic oscillations due to V® eT . 

We now proceed to obtain a priori bounds on "JP" 1 allowing to employ Theorem ll.4l An 
upper bound to the energy is easily derived by taking a trial function of the form 

trial CTm)* 



(4.4) 



\\UF n cr {-/m)y\\ L2{R3N) 

with ^ e H 1 (M. 3N ) normalized in L 2 (M. 3N ) independent of m. Using the results of Section[5] 



as above we obtain 



£ m [* trial ] = Nm- l E% T + C ol [*](l + O(m)). 



(4.5) 



As ^ does not depend on m we can use (|1.12j) to replace the complicated interaction F crys 
by the generalized Pekar interaction Ff m 



£ m [^] = iVm- 1 ^ + £f M [*](1 + 0{m)) + o(l). 

As this holds for any 'J g H 1 (M 3N ), we clearly have proved that 

lim sup E m (N) < NE? C * + EfjN). 
ra->0 



(4.6) 
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A rough lower bound to £ m [^ m ] is obtained by dropping the Coulomb term in (|4.2p and 
using Lemma 11.11 to bound the interaction with the crystal from below : 



juppcr / ™ i r i , >\ 

£ m [*m] > + ( E 2 J M3N l V ^™| 2 - 2° { m3 P^f(rn-),rn 3 p^(m-))\(l + 0(m)) 



NE? y ? 



1 \ J r3n iV^r'l 2 + \d (p^,p^fj (1 + 0(m)) 

where Ci and C2 are two positive constants. The third term in the right-hand side of (|4.7|) 
is given by the Pekar functional with a non optimal constant and is thus bounded below 
independently of to. Combining with the upper bound (|4.6|) we infer 



P^P°' ! P^P° l ) < C. 



(4.8) 



In the one-body case, this reduces to 



Iv^f + ^d^lM^ 01 ! 2 ) <c, 



i.e. is a bounded sequence in iJ^R 3 ). In the case N > 2, the Hoffmann-Ostenhof [21] 

inequality 



allows to deduce from (|4.8[) that (^/p $po i) is a bounded sequence in H 1 (R 3 ). 

We can now apply Theorem 11.41 to the sequence (^p? 1 ) when N = 1, or (^/p^poi) when 
iV > 2, and obtain from (O and (O 

£ m [* m ] = Nm- l E>? + C: [^ ol ](l + 0(m)) + o(l). 

On the other hand it is not difhcult to see that 



l^m 1 1 ^2 (jj3iv -j 1 + O(to) 



so, writing 



we have 



fl/Pol _ ||vl/P ol || Vj/Po! 

II m Hl 2 (R 3jv ) m 



f m [$ ra ] = JVm- 1 £r+£ i 



p 

EM 



pol 



(1 + 0{m)) + o(l). 



(4.9) 



Recalling f|4.6[) . we deduce that (^j^ ) is a minimizing sequence for Ef (N) and the conclu- 
sions of Theorem 11.31 follow by using Theorems 11.11 and 11.21 □ 
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